Entanglement of movable mirror and cavity field enhanced by an optical parametric amplifier
Zhang Cai-yun†, , Li Hu, Pan Gui-xia, Sheng Zong-qiang
Anhui University of Science and Technology, Huainan 232001, China

 

† Corresponding author. E-mail: zcylh9@163.com

Project supported by the National Natural Science Foundation of China (Grant No. 11247001), the Scientific Research Foundation of the Higher Education Institutions of Anhui Province, China (Grant No. KJ2012A083), and the Doctor (Master) Fund of Anhui University of Science and Technology, China.

Abstract
Abstract

A scheme to generate entanglement in a cavity optomechanical system filled with an optical parametric amplifier is proposed. With the help of the optical parametric amplifier, the stationary macroscopic entanglement between the movable mirror and the cavity field can be notably enhanced, and the entanglement increases when the parametric gain increases. Moreover, for a given parametric gain, the degree of entanglement of the cavity optomechanical system increases with increasing input laser power.

1. Introduction

Entanglement is one of the most attractive topics of quantum mechanics, which has been widely applied in diverse embranchments of physics. It is of basic and practical significance that we can generate entanglement among mesoscopic and even macroscopic systems. With the advance of technology, the Schrödinger cat state which describes the superposition states of macroscopic systems has been realized in experiments.[1] A macroscopic entanglement state can be prepared by optomechanical coupling via radiation pressure. Since the larger number of photons can lead to stronger radiation pressure, enhancing the photon number to increase the radiation pressure can acquire robust entanglement of optomechanical systems. A cavity optomechanical system is an available candidate system to study the entanglement of macroscopic states. Recently, there has been a growing interest in investigating entanglement in macroscopic systems.[213] In optomechanical systems, the quantum measurement theory and its applications have been explored,[3] a large variety of nonclassical states of both the cavity field and the movable mirror can be generated,[4] and the robust stationary entanglement can be produced between the intracavity mode and the mechanical mirror.[6] Reference [8] obtained a steady state entanglement of the motion of two dielectric membranes, which are suspended inside a Fabry–Perot cavity. The properties of optomechanical entanglement in a coupled cavity array with a movable mirror have been studied in Ref. [10]. Barzanjeh et al.[11] proposed a scheme for the realization of stationary continuous-variable entanglement in a hybrid tripartite system formed by an optical cavity and a microwave cavity, both interacting with a mechanical resonator. Chiara et al.[12] investigated the mirror–light entanglement in a hybrid optomechanical device formed by a Bose–Einstein condensate. In recent years, the experiments have realized the macroscopic entanglement.[14,15]

With the improvement of the qualities of the nonlinear crystals, optical parametric amplifiers (OPAs) including degenerate OPA and nondegenerate OPA, have shown tremendous applications in generating squeezed and entangled states.[1624] Agarwal[17] investigated interferences in the quantum fluctuations of the output of an OPA, and demonstrated that the interferences can be manipulated by choosing the squeezing of the input field. Reference [23] obtained broadband entangled light through cascading nondegenerate OPAs, and the application of the entangled light from the cascading nondegenerate OPAs to broadband teleportation has been discussed.

Recently, the enhancement of the entanglement has been investigated by introducing nondegenerate OPAs.[2527] Chen et al.[25] theoretically showed that the degree of the correlation of the input entangled beams can be improved by the nondegenerate OPA inside an optical cavity. In Refs. [26] and [27], it was experimentally demonstrated that the entanglement degree of the entangled state can be improved and manipulated by using nondegenerate OPAs. In this paper, we propose a scheme to produce the stationary macroscopic entanglement between a cavity field and a movable mirror in an optomechanical system filled with OPA. With the existence of the OPA, the entanglement can be pronouncedly increased compared to that without OPA, and the entanglement increases with the increasing parametric gain of the OPA. Furthermore, for a fixed parametric gain, the larger entanglement of the optomechanical system corresponds to the larger input laser power.

2. Model and steady-state entanglement for the system

As shown in Fig. 1, we consider an OPA confined inside an optical cavity with a fixed mirror and a perfectly reflective movable mirror driven by a laser. The Hamiltonian of the system can be given as

The first term describes the energy of the cavity, where a (a) is the annihilation (creation) operator of the cavity fields with frequency ωc = 2πc/L, and L is the cavity length. The second term is the energy of the mechanical oscillator, where Q (P) is the dimensionless position (momentum) operator of the mechanical oscillators with frequency ωm. The third term represents the radiation pressure with the coupling rate

The fourth term gives the coupling between the input laser field and the cavity mode, where ɛ is related to the input power by and ωL is the frequency of the driving laser. The last term is the coupling between the OPA and the cavity mode, where G is the nonlinear gain of the OPA and θ is the phase of the field driving the OPA. In the interaction picture, the Hamiltonian is given as

Fig. 1. Sketch of the system. An optical parametric amplifier is inserted in a Fabry–Perot cavity with a fixed mirror and a movable mirror.

Considering all of the noise and damping terms, we can investigate the dynamics of the system by using the quantum Langevin equation. By introducing phenomenologically the damping and noise terms into the Heisenberg operator equation, the quantum Heisenberg–Langevin equation can be obtained as

where κA is the decay of the cavity field and Ain describes the input noise. From Eq. (2), we have

where Δ = ωcωL +χQ. The quantum Brownian noise operator ξ with zero mean value comes from the coupling of the movable mirror to the environment, and has the following correlation function

where kB is the Boltzmann constant, and T is the temperature of the mechanical mirror. When ωm/γm ≫ 1, the correlation function of the noise ξ(t) can be written as 〈ξ(t)ξ(t′)+ξ(t′)ξ(t)〉/2 ≅ γm(2N + 1)δ (tt′), where N = [exp(ħωm/kBT) −1]−1 is the mean thermal excitation number of the resonator. The noise operator ain satisfies 〈ain(t)〉 = 0 and the nonzero time-domain correlation functions , , where n = [exp(ħω0/kBT)−1]−1 is the equilibrium mean thermal photon number.

In order to study the steady state entanglement of the system, we adopt the standard methods of quantum optics[28] to solve Eq. (4), namely, the steady state solutions of the system can be obtained by setting all the time derivatives of Eq. (4) to zero. Consequently, we obtain the steady state values of the system

When the cavity is very intensely driven, i.e., as ≫ 1, each operator of the system can be expanded as the sum of its steady-state mean value and a small fluctuation with zero mean value, namely, the operator can be written as ô = os+δô, and the dynamics of the small fluctuations is around the steady state of the system. We have linearized Langevin equations as follows:

where we have defined quadrature operators for the cavity field as

Equation (6) can be written in a compact form

where

and

are the column vectors of the fluctuations and the noise sources, respectively. The drift matrix A is written as

When the real parts of all the eigenvalues of matrix A are negative,[6,29] the system is stable and reaches its steady state as t → ∞. It is not easy to obtain analytically the eigenvalues of the drift matrix, so we will guarantee steady state conditions via a numerical method. We define the element of the covariance matrix of the quantum fluctuations Vij(∞) = [〈fi(∞)fj(∞) + fj (∞)fi(∞)〉]/2, and investigate the nature of the steady state of the cavity optomechanical system. When the system is stable, the steady-state correlation matrix satisfies the Lyapunov equation[30]

where D = diag[0, γm(2n + 1),κ (2N + 1),κ (2N + 1),0,0] is the noise correlation matrix. We can directly solve the covariance matrix.

We aim to investigate the properties of entanglement of the cavity optomechanical system containing OPA. We employ the logarithmic negativity to quantify the entanglement. The logarithmic negativity can be defined as[31]

where with Σ(V) = det A + det B − 2det C, and it is the smallest symplectic eigenvalue of the partially transposed covariance matrix given by

According to the definition of Eq. (11), the system is entangled if EN > 0, or equivalently, η < 1/2.

We now discuss the properties of the entanglement of the cavity optomechanical system. The entanglement of the system is displayed in Figs. 2 and 3, and the logarithmic negativity EN is plotted as a function of the normalized detuning Δ/ωm. In Fig. 2, we show the logarithmic negativity of the system for three values of G, i.e., G = 0, 0.8κ, and 0.9κ. From Fig. 2, we notice that the larger entanglement between the movable mirror and the cavity mode can be obtained when the OPA exists, and the entanglement increases with the increase of G. This is because increasing the parametric gain makes the coupling between the cavity field and the movable mirror stronger due to the increase of the photon number inside the cavity. The photon number inside the cavity is 6.29 × 108, 8.02 × 109, and 10.34 × 109 for G = 0, G = 0.8κ, and G = 0.9κ, respectively. Moreover, in Fig. 2, it can be observed that the entanglement with the OPA (dashed and solid lines) can be notably enhanced compared to that without OPA (dashed-dotted line) between the movable mirror and the cavity field.

Fig. 2. Entanglement between the cavity field and the mechanical mirror, and the dashed-dotted, dashed, and solid lines correspond to G = 0, 0.8κ, and 0.9κ, respectively. For all plots, the parameters are L = 5 mm, m = 20 ng, ωm = 2π × 107 Hz. The laser has wavelength λ = 810 nm and power = 10 mW, the mechanical quality factor Q′ = ωm/γm = 6700.

Figure 3 describes EN of the cavity optomechanical system for different input laser power at G = 0.8κ. From Fig. 3, it can be observed that the entanglement of the optomechanical system increases with the increase of the laser power when parametric gain G is fixed. The reason is that the coupling between the cavity field and the mirror is increased with increasing input laser power due to the increase of the photon number.

Fig. 3. Entanglement between the cavity field and the mechanical mirror when G = 0.8κ, the dashed-dotted, dashed, and solid lines correspond to = 7 mW, 10 mW, and 15 mW, respectively. The other parameters are the same as those in Fig. 2.
3. Discussion and conclusion

We give a brief discussion on the experimental feasibility of our program. To numerically calculate the logarithmic negativity of the optomechanical system, we have chosen feasible parameters based on the present experiment.[32] In fact, it is an experimental challenge to detect the entanglement of macroscopic mechanical systems. However, some promising schemes have been proposed[6,33] to relatively easily detect the quantum correlation by homodyne measurement techniques, so quantum entanglement can be indirectly detected by quantum correlation. We believe that our scheme is realizable in experiments with the development of science and technology.

In conclusion, we have proposed a scheme to generate stationary macroscopic entanglement of an optomechanical device filled with OPA. Our study shows that the entanglement of the cavity optomechanical system can be pronouncedly enhanced due to the introduction of the OPA, and the larger entanglement corresponds to the larger parametric gain G. Moreover, for a given G, the degree of entanglement increases when the input laser power increases.

Reference
1Brune MHagley EDreyer JMaître XMaali AWunderlich CRaimond J MHaroche S 1996 Phys. Rev. Lett. 77 4887
2Law C K 1995 Phys. Rev. 51 2537
3Chen XLiu X WZhang K YYuan C HZhang W P 2015 Acta.Phys. Sin. 64 164211 (in Chinese)
4Bose SJacobs KKnight P L 1997 Phys. Rev. 56 4175
5Bose SJacobs KKnight P L 1999 Phys. Rev. 59 3204
6vitali DGigan SFerreira ABöhm H RTombesi PGuerreiro AVedral VZeilinger AAspelmeyer M 2007 Phys. Rev. Lett. 98 030405
7Mancini SGiovannetti VVitali DTombesi P 2002 Phys. Rev. Lett. 88 120401
8Hartmann M JPlenio M B 2008 Phys. Rev. Lett. 101 200503
9Paternostro MVitali DGigan SKim M SBrukner CEisert JAspelmeyer M 2007 Phys. Rev. Lett. 99 250401
10Wu QXiao YZhang Z M 2015 Chin. Phys. 24 104208
11Barzanjeh SVitali DTombesi PMilburn G 2011 Phys. Rev. 84 042342
12Chiara G DPaternostro MPalma G M 2011 Phys. Rev. 83 052324
13Tan H TLi G X 2011 Phys. Rev. 84 024301
14Thompson J DZwickl B MJayich A MMarquardt FGirvin S MHarris J G E 2008 Nature 452 72
15Verhagen EDeléglise SWeis SSchliesser AKippenberg T J 2012 Nature 482 10787
16Collett M JGardiner C W 1984 Phys. Rev. 30 1386
17Agarwal G S 2006 Phys. Rev. Lett. 97 023601
18Mehmet MVahlbruch HLastzka NDanzmann KSchnabel R 2010 Phys. Rev. 81 013814
19Eckstein AChrist AMosley P JSilberhorn C 2011 Phys. Rev. Lett. 106 013603
20Zhang JYe CGao FXiao M 2008 Phys. Rev. Lett. 101 233602
21Zhao C Y 2015 Chin. Phys. 24 040302
22Silberhorn CLam P KWeiss FKönig OKorolkova NLeuchs G 2001 Phys. Rev. Lett. 86 4267
23He W PLi F L 2007 Phys. Rev. 76 012328
24Yan Z HJia X JSu X LDuan Z YXie C DPeng K C 2012 Phys. Rev. 85 040305
25Chen H XZhang J 2009 Phys. Rev. 79 063826
26Shang Y NJia X JShen Y MXie C DPeng K C 2010 Opt. Lett. 35 853
27Zhou Y YJia X JLi FYu JXie C DPeng K C 2015 Scientific Reports 5 11132
28Walls D FMilburn G J1998OpticsBerlinSpringer10.1007/978-3-540-28574-8
29Dejesus E XKaufman C1998Phys. Rev. A355288
30Genes CMari ATombesi PVitali D 2008 Phys. Rev. 78 032316
31Adesso GSerafini AIlluminati F 2004 Phys. Rev. 70 022318
32Marquardt FChen J PClerk A AGirvin S M 2007 Phys. Rev. Lett. 99 093902
33Mazzola LPaternostro M 2011 Phys. Rev. 83 062335